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Abstract 

In the limit of small values of the aspect ratio parameter (or wave 
steepness) which measures the amplitude of a surface wave in units of its 
wave-length, a model equation is derived from the Euler system in infi- 
nite depth (deep water) without potential flow assumption. The resulting 
equation is shown to sustain periodic waves which on the one side tend 
to the proper linear limit at small amplitudes, on the other side possess 
a threshold amplitude where wave crest peaking is achieved. An explicit 
expression of the crest angle at wave breaking is found in terms of the 
wave velocity. By numerical simulations, stable soliton-like solutions (ex- 
periencing elastic interactions) propagate in a given velocities range on 
the edge of which they tend to the peakon solution. 
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1 Introduction 

The description of the propagation of surface waves in an ideal incompressible 
fluid is still a classical subjet of investigation in mathematical physics as no 
definite comprehensive answer to the problem has been given yet. In the limit 
of shallow water, surface gravity waves have been intensively studied and many 
model equations were introduced by various approaches, with great success. The 
nonlinear deep water case is more cumbersome and there does not exist today a 
simple model as universal as the shallow water equations (Korteweg-de Vries or 
Boussinesq) which would result from an asymptotic limit of the Euler system. 

The inherent technical differences between shallow and deep water are mainly 
due to the fact that the two natural small parameters used for perturbative 
analysis of the Euler system in shallow water loose their sense in the deep water 
case (depth h — > oo). Indeed, these parameter are ei = a/h, which measures 
the amplitude a of the perturbation scaled to the fluid depth h, and £2 = ^^/A'^, 
which measures the depth in units of wavelength A. 
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By perturbative expansion in ei, and fixed finite 62, one obtains the nonlinear 
shallow water equation Retaining ei and €2 (but not their product) leads 
to different versions of the Boussinesq equations P|, from which the Korteweg- 
de Vries equation is derived by assuming a small amplitude wave moving in a 
given direction All these model equations govern asymptotic dynamics of 
long wavelength wave profiles. 

One way to obtain a small parameter in deep water is to take into account 
that deep water waves typically result from a superposition of wave components 
more or less close to a fundamental carrier wave. The small parameter then 
measures the variations of the envelope scaled to the one for the carrier wave. 
Such nonlinear modulation of wave trains is worked out perturbatively by means 
of slowly varying envelope approximation (SVEA) which usually leads in 1+1 
dimensions to the nonlinear Schrodinger model ^ Ej. For a full account on 
modulation of short wave trains on water of intermediate or great depth we refer 
to ini E] . The procedure provides nonlinear dynamics of surface waves under 
the form of nonlinear modulation, and the drawback is that the dynamics of the 
wave profile itself remains unknown. 

Based on the theory of analytical functions and perturbation theory, a model 
for the profile of the free surface wave in water of finite depth, involving the 
Hilbert transform operator, was derived in 0. Although this model possess 
a well-defined deep water limit, the resulting equation cannot be studied by 
known techniques to compare it to KdV-like models. Other model equations, 
built to fit the properties of waves on deep water can be found in EI| . Their 
dispersion relations coincide exactly with that of the water waves on infinite 
depth but their nonlinear terms are chosen ad hoc to reproduce Stokes waves. 

Our purpose here is to study the asymptotic dynamics of the very profile 
of a surface wave in deep water in the weak nonlinear limit by assuming a 
dependence on the vertical coordinate close to the linear harmonic one. The 
dispersion relation of a fluid on a depth h 

LO^ ^ gktanh{kh), (1.1) 

leads for long waves on shallow water (parameter kh small) to the nondispersive 
relation 

uj = ky/gh Vp = Vg = \/gh, (1-2) 

where Vp is the phase velocity and Vg the group velocity. However, waves on 
deep water (parameter kh — > cjo) are dispersive as from Hl.l|l 

u; = y^fc ^ Vp = Vg = (1.3) 

This deep water dispersion relation will be one of our main guide in the process 
of finding a limit model whose linear limit is constrained by p.3|l . 

Our approach follows the method of Green and Naghdi ^2 ^] for surface 
waves in shallow water which assumes an anzatz for the dependence of the ve- 
locity components on the vertical dimension z. This anzatz does not produce 
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an exact solution of the full Euler system and the game consists in replacing 
one of the equations with its integrated expression. This comes actually to mak- 
ing an average over the depth, which can be performed with different weights. 
Although weight is not determinant in the shallow water case we shall see 
that its choice is prescibed by consistency requirement at the linear limit. 

A limit model is then obtained by defining a small parameter which measures 
the amplitude of a surface wave in units of its wave-length, we call it the aspect 
ratio parameter (ARP). The approach combines the asymptotic analysis a la 
Whitham with the already standard method of multiple scales |14l 1151 IT^ I17| . 
This will be shown to lead to the following model 

13 5 1 

m - Vxxt - -jpxxx + -j^lx + VVx = ^VxVxx + -^VVxxx- (1-4) 

for the dimensionless deformation ri{x, t) of the free surface of deep water. 

The paper is organized as follows. In section [5] we introduce the Euler 
equations, their nondimcnsional version and the anzatz which, together with 
convenient average, enables to reduce the initial three-dimensional problem to 
a two-dimensional one. The resulting model, equations (|2.25f) below, is in fact 
the deep water version of the celebrated Green-Naghdi system for surface waves 
in shallow water [111 [T^ . 

In section^ the nonlinear dispersive model (|1.4|) that governs propagation of 
waves of small ARP in deep water is derived by the method of multiple scales and 
perturbative expansion. Expressed in physical units, the model appears with 
fc-dependent coefficients, reminiscent of what occurs in the SVEA approach of 
deep water modulation. 

In section 0] analytical and numerical analysis of the progressive periodic 
waves are performed. Numerical simulations show the existence of periodic 
waves than tend to peak as amplitude grows. The crest angle at the peak is 
calculated and we obtain an explicit expression in terms of the wave velocity. 
Moreover numerical simulations show also the existence of soliton-like solutions 
(though we cannot provide analytic expression) that tend, at the edge of the 
allowed velocity range, to the peakon solution (who is given an explicit analytical 
expression) . 

2 The model equations in deep water 
2.1 General settings. 

Let us consider the Euler equations in physical dimensions where the particles 
of the fluid are identified in a fixed rectangular Cartesian system of center O and 
axes {x,y,z) with Oz the upward vertical direction. We assume translational 
symmetry in y and thus consider a sheet of fluid in the xz plane. The velocity 
of the fluid in this plane is the vector {U, W) where U{x, z, t) is the horizontal 
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component and W{x,z,t) the vertical one. This fluid sheet is moving on a 
bottom at z = — oo and its free surface is z = r]{x, t). 

The continuity equations and the Newton equations in the flow domain read 

U, + W,= 0, (2.1) 

pU = -Px, (2.2) 

pW = -P, ~ gp, (2.3) 

where P{x,z,t) denotes the pressure, p the uniform density and g is the grav- 
ity. Subindices mean partial derivatives and overdot means material derivative 
defined as usual by 

F = Ft + UF, + WF,. (2.4) 

The boundary conditions at z = — cxd simply state that the velocity component 
vanish (U and — > 0), while boundary conditions at the free surface 
z = ri{x, t) state that P is the atmospheric pressure Pa and that the total 
derivative of the surface equation z — rj = Q vanishes, namely 

z^-q : rjt + Ur]^ - W = 0. (2.5) 

We are interested here in finding the evolution equations for the free surface 
by studying the nonlinear deformation of a particular linear wave profile with 
given arbitrary wave number k. 



2.2 Dimensionless Euler system. 

The linear progressive monochromatic solution of the linear limit of the above 
Euler system reads fTHI 

C^iinoar = C^o cos{kx - ujt) e'"^, oj"^ = gk, (2.6) 

where the dispersion relation is indeed of the deep-water class (|1.3|l . 

Thus, given the wave number fc, we may scale the original space variables x 
and z with k, the time variable with \/fcg, the velocity components U and W 
with \/k/g, and the pressure P with k/pg. For convenience we keep the same 
notations for the adimensional variables and the Euler equations for z € [— c», 77] 
then become 

U^ + W,^ 0, (2.7) 
U = -P., (2.8) 
W^~P,~ 1, (2.9) 

whith the boundary conditions 

0, (2.10) 
(2.11) 

C/??., (2.12) 



-> —00 : 


U 


= w 


z = rj : 


P 


= Po, 


z — rj : 


Vt 


= w 
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where Pq = Pak/pg. Note that, as z is scaled with k, the dimensionless surface 
profile 77 is also scaled with k. 



2.3 Vertical velocity profile. 

Inspired thus by we assume an exponential vertical dependence of the 

horizontal velocity component and seek an approximate solution of the Euler 
system by starting with the anzatz 

U{x,z,t) = u{x,t) e\ (2.13) 

It does satisfy the boundary condition (|2.10() and gives 

W{x,z,t)^~u,{x,t) e', (2.14) 

obtained from the continuity equation H2.7|l and the boundary condition (|2.10(l . 
By computing the total derivative of the above velocity components we get 

[/ = ^ite^ (2.15) 
W = -u^^t e"" - uu^^ e^"" + ul e^"" , (2.16) 

which allows to calculate explicitely the pressure P out of (|2.9|l and boundary 
condition H2.11|l as 

P - Po = (r; - z) - u^t{e^ - e') + ^(u^ _ uu^,){e^^ - e^^). (2.17) 

This pressure diverges as z ^ —00 in Archimedean way as it must. 
The boundary condition (|2.12l) finally provides the evolution of 77 as 

+ {u e"), = 0. (2.18) 

To that point the anzatzs H2.13|l and H2.14|l , the expression 1)2. 17|) of P and the 
evolution H2.18|l of the free surface 77 satisfy exactly the differential equations 
(|2.7|l and H2.9|) together with the whole boundary conditions (|2.1U|) , H2.11|l and 
H2.12|l . The remaining equation to take into account is then the Newton's law 
H2.8|l which is of course not satisfied globally by the above expressions of W and 
P. 



2.4 Averaging Newton's law. 

The Newton's law is now taken into account through an average on the full 
depth with the weight e"^ which regularizes the diverging Archimedean term in 
the expression of the pressure (|2.17|) . This is explicitely 

dz e°'\U + P^)^Q. (2.19) 
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which eventually furnishes with use of (|2.15|) and H2.17|l 



, wt = -7^ ^ — ri-^^i + 7 — r^["^:r^ - "J] ■ (2-20) 

a + l a ox \ a + \ [a + 2) a J 

It is a model- dependent system where the solutions now depend on the external 
parameter a. 

The value of parameter a is fixed by demanding that the linear limit of the 
system (I2.18f)(l2.20f) possess solution with phase and group velocities given by 
p. 3(1 . In dimensionless units, these velocities are then required to be 

Vp^l, vg^]^. (2.21) 

In the linear limit of system H2.18|l (|2.2(J|) , rj can easily be eliminated to arrive at 

autt - (a + l)uj,x - Uxxtt = 0, (2.22) 
which possess the solution 

u = uq cos{qx — ojt) , uj'^{q) — — . (2.23) 

a + q"^ 

Requiring relation H2.21|l for the above dispersion law Lo{q) comes to require the 
relations 



uj{q) 



duj 



= \- (2.24) 



The first of these equations is verified for any a while the second holds if and 
only if a = 1 . 

As a result the basic system of equation is (|2.18|l and (|2.20|l with a = 1 
namely 

\ut e^" + dx (^e" - Uixt e^" + ^e^^[ul - uu^S) = 0, 

m + {u e"), = 0. (2.25) 

The system (|2.25|l for the couple of variables u{x^ t) and rj{x, t) is the net result 
of assuming the z-dependence (|2.13|l in the Euler equations and of making a 
convenient average of Newton's law on the depth (— oo, rj\. As displayed above, 
the crucial step is averaging H2.8|l with the precise weight e^. Indeed, while it 
can be proved from |13j that weighting the average is not determinant in the 
shallow water case, weight is essential here and the choice is justified by 
requirement ()1.3|l . 
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3 Small aspect ratio waves in deep water 



3.1 Generalities. 

The small parameter of the perturbative analysis of the system H2.25|l is, in 
dimensionless units, the maximum amplitude e of the free surface deformation 
ri{x,t). So we define 

T]{x,t) ^ eH{x,t), (3.1) 

and it is instructive to understand the physical meaning of e by returning to the 
physical units for which 

e = afc, (3.2) 

where a is the wave profile maximum amplitude and k the chosen wave number. 
Thus e is the parameter that measures the ratio of wave height to wave length, 
we call it the aspect ratio parameter (ARP). Then a small ARP (e < 1) means 
a flat deformation of the surface without reference to the absolute amplitude or 
to the depth. 

The dimensionless linear dispersion relation (w^ = 1) undergoes deviations 
due to nonlinearity (Stokes' hypothesis) which can be taken into account through 
the expansion 

Lo^ = l + e + e^ + ■■■ . (3.3) 
Therefore the phase expands as 

X - Lut ^ X - t - et - e'^t -\ , (3.4) 

which suggests the following definition of the new variables y, r, • • • , 

y = X — t, T — et, V = e^t, ■ ■ ■ . (3.5) 

The function H{x,t) is then represented as the function H{y,T,h', ...,) of the 
new variable y and the slow variables r, v, • • • , for which hold the derivation 
rules 

dx ^dy, (3.6) 
dt^- dy + tdr + + • • • . (3.7) 

These rules are now explicitely applied step by step to our system ()2.25|) . 

3.2 Asymptotic expansion. 

Since exp(eiJ) = l + eiJ + 0(e^) we can consider system (|2.25f) at order e°, next 
at orders e'^ and e and so on. In terms of H , y and r we get 

i-dy + tdr + Oie^))eH + [u(l + eH + Oie^))]y = (3.8) 
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This equation gives u in term of derivatives and integrations of H 

u^eH -e^{H^ + T Hrdy')+0{e^) (3.9) 

J —OQ 

where it is not assumed that u ov H (and derivatives) vanish as ?; — s- — oo. 
Instead the above equation readily gives the usefull relation 

u(-cx),t) ^ eH{-oo,T)[l - eH{-oo,T) + 0{t^)]. (3.10) 

Now equation H2.18|l is more convcneintly written in the form 

\a+^B-C-]^D = (3.11) 

with A,B,C and D given by 

A={u,te^'^)^, S = ((mm,, - (u,)2)e3'7)^, C = (e"),, D ^ ute''\ (3.12) 
Using the operator expansion 

dl,^~dly + edl, + Oie^) (3.13) 

we obtain 

A^-^{{l + 2eH + 0{e^))x 
ay 

J — OO 

which at orders e and provides 

A = -eHyyy + e^AHyHyy + 2Hyyr) + O^e^). (3.14) 

Coming now to the computation of B wc may write 

B=-^{[\ + ieH + 0{e^)]x 

oy 

[(ei/ _ e2ij2 _ ,2 r ^^^y, ^ oie')){eHyy + ©(e^)) _ e^Hyf + ©(e^)]} 

J —OO 

which at order provides 

B = e^HHyyy - HyHyy) + O (3.15) 

Finally C and D at order are given by 

C He'")y = ^Hy + e^HHy + ©(e^), (3.16) 
D^[~dy + edr+0(e^)][eH-e'{H^+ T i/,dy') + 0(e^)] x 

^ — OO 

[l + 2ei/ + 0(e2)] 
= -eHy + e^2Hr+0{e^). (3.17) 
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We substitute A, B, C and D in H3.11|l . keep the terms in e° and to obtain 
eventually 

10 2 

-Hy - Hyyy + e{—HyHyy + -HHyyy - 2HHy + 2Hyyr ^ 2H r) = (3.18) 

or in the original (dimensionless) x and t variables 

In terms of the free surface rj = eH the above equation reads 

13 5 1 ,^ 

m - Vxxt - -jTj^xx + -jVx + VVx = -^VxVxx + -^llVxxx, (3.20) 

as announced in the introduction. This is the equation which describes asymp- 
totic nonlinear and dispersive evolution of small ARP waves in deep water. 



3.3 Discussion. 

We note first that the evolution l|3.20|l of the free surface, written in physical 
dimensions, namely by 77 — > fcry, x kx and t — > t^/gk, reads 

vt-^ Vxxt -^sji Vxxx+^J ^vx + ^/gk mx 



^ yf '■^''^^^^ ^ VVxxx)- (3.21) 

It is therefore a fc-dependent equation that describes the nonlinear deformations 
of the wave of given k. 

Next, by defining 

F(x,<) = iry(x,i) + i (3.22) 



we transform the equation (|3.2U|) into 

Ft Fxxt "i" '^FFx — ^FxFxx ^~ FFxxx- (3.23) 

Equations of the form 

Ft ~ Fxxt + {b+ l)FFx = bFxFxx + FFxxx, (3.24) 

are completely integrable only in the cases 6 = 2 (Camassa-Holm equation) [2D] 
and 6 = 3 (Degasperis-Procesi equation) j21[l22[l2n] . Our equation H3.23|l does 
not enter this class and we cannot tell more about integrability. 

Finally, from the above perturbative expansion we also have, according to 
H3.9|l the same equation for the variable u{x,t), which can also be written as 
the conservation law 

1 3 12 1 

dt{u - Uxx) = dx{-Uxx " 2" ~ ^ 3^^ ^ -^uuxx)- (3.25) 
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4 Progressive waves 



4.1 Periodic progressive waves. 

The waves of constant profile propagating at velocity c are the solutions of H^12()|l 
for 

Vix,t)=r^{0, ^ = x-ct, (4.1) 
which implies the ordinary nonlinear differential equation 

3 1 5 1 

(2 ~ ^ ^2 ~ '^''^'^^^ ^ " 3^?^« + 3 (4.2) 

which can be integrated once to give 

3 1 12 1 

(2 - c)?7 - (2 - + 2^^ 3^1 ^ '^^^ii ^' 

where A is an integration constant to be determined from the chosen boundary 
conditions. 

To study the periodic solutions, we integrate (|4.2|l numerically under the 
initial values 

77(0)^0, rye(0) = 0, r/55(0) = «, (4.4) 

and obtain a typical set of profiles displayed in figure^for different values of the 
parameter a. Such a simulation shows the wave crest peaking, together with 
wavelength variations, as wave amplitude increases. We have checked that low 
amplitude profiles (obtained for small values of a) do tend to the linear solution 
a(cos(0 + !)■ 




Figure 1: Profiles of the dimensionless surface deformation rj solution of (|4.2|l 
with (1131) for c = 1 and a taking the values 0.2, 0.3, 0.4, 0.446, 0.4499 and 0.45, 
from bottom to top. 

The value of the maximum amplitude rjs reached at wave peaking can be 
infered from the equation 1)4. 3|l . it is the value for which the coefficient of the 
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second derivative vanishes, namely 

Vs=5{c~^). (4.5) 

We have found actually that this occurs for the initial data (|4.4() with the thresh- 
old Us given by 

«. = |j(^-c). (4.6) 
At the above threshold values, the crest angle can be evaluated as follows. 



4.2 Crest angle. 

First the integration constant A of l|4.3|l is evaluated at ^ = by means of (|4.4|l 
and H4.fi|l . we obtain 

^->-^)(^c). (4.7) 

Next we evaluate equation H4.2|l at the limit ^ — > (limit to the left ^7 or to 
the right ^+), where is the value for which the threshold is reached, namely 
vi^s) = Tjs- We obtain this way the limit of the second derivative by 

6 

liin (????) = 7"- (4-8) 




Figure 2: Values of the crest angle, in degrees, at wave peaking for velocities c 
in the range [0.6, 1.2]. 



To get the limits of the first derivative we evaluate now equation (|4.3|l at 
^ = ■ With help of the above value of A we obtain 

^lhn(77|) = ^(c-i)(^c+l), (4.9) 
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and thus the limit of 77^ in can be either positive or negative. It is then clear 
that a bounded periodic solution like the one displayed on figure ^ implies the 
following solution (remember c > 1/2) 



limjry^) = - lim (7?^) = 



1/2 



(4.10) 



Last, from the above value of the slope at ^s, the crest angle at wave peaking 
reads 



9 = 2 arctan • 



1 6 

2) ^r- 



-1/2' 



(4.11) 



and the graph of this expression is plotted on the figure El At c = 1 the angle 
is 6* = 65 d°. 

It is worth remarking that the analytic expression of the crest limit angle, 
which we have cheked on numerical simulations, actually reults in a part from 
the phenomenological expression H4.6|l obtained by observations of numerical 
simulations of equation l|4.2|l with boimdary data H4.4(l . 



4.3 Peakon solution. 

The coefficients of the ODE (|4.2f) show that periodic solutions hold if c > 1/2 
and a threshold is reached at c = 5/2 for which waves of vanishing amplitudes 
are already singular. At this value, the model supports the peakon solution 

77(e)-6exp[-|^|/N/2], e = (4.12) 

Remarkably enough, this peakon also appears as the limit of the solitary wave 
solution when its velocity reaches the value 5/2. 



4.4 Soliton-like waves. 

Numerical simulations has revealed the existence of stable soliton-like excita- 
tions that travel without deformation at constant speed in the range [3/2, 5/2]. 
These are performed by observing that 



is an approximate solitary wave solution for (|1.4() according to 



4'' - Vil - Ivit + Ivi'^WM'' - Iv^^V^ - Iv'-'vl^. = OinSl). (4.14) 



Hence defining 



3 S 

- + 5 ^ {7jo = 3S, 4^2^^}, (4.15) 
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expression (|4.13|l tends to the exact solution for S ^ (for which the amphtude 
770 vanishes). 

We have performed a series of numerical simulations of (|1.4|) by first writing 
it in the frame moving at velocity v, namely 

r?t + - v)ilx - Vxxt + (u - ^)rixxx + VVx = ^VxVxx + ^mxxx- (4-16) 
and then by using the following initial-boundary value problem on a; G [0, i] 

rKO,t) =r,,(x,0), 77(0,0=0, r7(£,t) = 0, 7fe(L,i) = 0, (4.17) 
with an initial position = 2L/3 and a velocity 

3 



+ 5, 5e]0,i[. 



(4.18) 



The method uses a standard second order implicit finite difference scheme with 
a mesh grid dimension 0.1. 




v=1.6 t=80 








50 100 150 



50 100 150 




v=1.9 t=300 



Figure 3: Solitary wave profiles in two examples with (5 = 0.1 (up) and 5 = 0.4 
(down) in the rest frame of the approximate solution H4.13|l . 

A few typical examples are displayed in fig|3| where we compare the solitary 
wave to the initial condition. As expected, when 5 « 1 (i.e. v ^ 1.5) the 
solitary wave resembles much the initial condition while for larger 5, a pulse 
reshaping occurs by means of emission of waves (at phase velocity 1). After 
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reshaping, the solution is remarkably stable. When S ^ 1 (i.e. v 5/2), the 
solitary wave tends to the peakon solution (|4.12() . 

The found solitary wave solutions behave like solitons do in an integrable 
system. Indeed, after reshaping, a two-soliton interaction effects only the po- 
sition of the pulses and does not alter the individual shapes. An example is 
diplayed in figure 01 




(SO so lOO 120 140 ftO SO 1 OO 1 20 1 40 60 SO 1 OO 120 140 



Figure 4: Solitary wave interaction in the rest frame of the soliton at velocity 
1.8, the other soliton having velocity 2. 



Conclusion and comments 

The asumed exponential vertical velocity profile has led to a simple model for 
surface waves on deep water. The derivation is rigorously made by perturbative 
analysis assuming waves with small aspect ratio parameter, which constitutes a 
novel approach of the problem. 

The resulting model H3.2Q(I sustains periodic waves which tend to peak as 
their amplitude increase and reach a threshold amplitude rjs given by (|4.5|l . At 
this value, the crest angle can be explicitely computed in terms of the velocity 
by H4.11|l . However, in the region of the peak, higher order derivatives diverge 
and the perturbative expansion does not hold anymore. Hence the precise value 
of the crest angle should be considered with care and the peaking only gives an 
indication of a behavior due to nonlinearity. 

As mentionned in the introduction, it is necessary to release the potential 
flow approach which is not consistent with the anzatz H2.1^{|l . Indeed, assuming 
the existence of a potential flow (j) by 

([/, ly) = V(/), (4.19) 
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the expressions (|2.13|) and (|2.14|) . together with (|2.7I) . readily provide 

(j){x,z,t) ^ ip{x,t)e'', Lpa,x + V = 0. (4.20) 

Although one may obtain nontrivial time-dependence from the Euler system in 
this particular context, the above linear x-dependence of ip{x,t) cannot catch 
nonlinear deformations in a consistent way. 

Last but not least, it is remarkable that equation H3.21|l possess /c-dependent 
coefficients, where k is the wave number of the selected particular wave. It is 
something new for a model describing a wave profile. Indeed asymptotic non- 
linear and dispersive long scales models for surface water waves (as Korteweg- 
deVries, modified KdV, Benjamin-Bona-Mahoney-Peregrine, Camassa-Holm) 
have fc-independent coefficients. Conversely, fc-dependent coefficients are com- 
mon in the context of modulation of wave trains, as in nonlinear Schrodinger, 
modified NLS, Davey-Stewartson, etc... 
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